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ABSTRACT 

We investigate Venus Express (VEX) observations of magnetic field fluctu¬ 
ations performed systematically in the solar wind at 0.72 Astronomical Units 
(AU), between 2007 and 2009, during the deep minimum of the solar cycle 24. 
The Power Spectral Densities (PSD) of the magnetic held components have been 
computed for the time intervals that satisfy data integrity criteria and have been 
grouped according to the type of wind, fast and slow dehned for speeds larger 
and respectively smaller than 450 km/s. The PSDs show higher levels of power 
for the fast than for the slow wind. The spectral slopes estimated for all PSDs 
in the frequency range 0.005-0.1 Hz exhibit a normal distribution. The average 
value of the trace of the spectral matrix is —1.60 for fast solar wind and —1.65 
for slow wind. Compared to the corresponding average slopes at 1 AU, the PSDs 
are shallower at 0.72 AU for slow wind conditions suggesting a steepening of the 
solar wind spectra between Venus and Earth. No signihcant time variation trend 
is observed for the spectral behavior of both slow and fast wind. 

Subject headings: interplanetary medium — magnetic fields — plasmas — solar wind 


turbulence 



3 


1. Introduction 


The spectral properties of the magnetic held and plasma huctuations in the solar wind 
have been investigated in-situ over several decades, for a broad range of frequencies and for 
various radial distances. It was found that the power spectral density (PSD) of magnetic 
held huctuations exhibits three diherent power law regimes, P{k) = characterized 


by diherent exponents: (i) a —1 for smaller k (e.g. Matthaeus & Goldstein, 1986, Verdini 


et al. 

2012), (ii) a ~ —5/3 for the intermediate k (e.g.. 

Marsch & Tu 

1990 


k is also anisotropic, the huctuations parallel and perpendicular to the magnetic held may 


exhibit diherent power law index (see, e.g., Dasso et ah, 2005 Horbury et ah, 2012) (hi) 


a < —2.5 with a minimum index close to —4.5 (Leamon et al. 1999 Bruno et al 2014b) for 


the largest k (see, also, Coleman, [1968; Stawicki et ah, 2001, and Bruno & Carbone, 2013 


Alexandrova et ah, 2013, for a review). [Frisch (1995) described the three characteristic 


power law regimes, separated by spectral breaks, as the magnetohydrodynamic equivalents 
of the scale ranges of the classical hydrodynamic turbulence: (i) the driving (or energy 
containing) range, (ii) the inertial range, dominated by nonlinear turbulent interactions 
that transfer the energy over multi-scales and (hi) the dissipation range. The physical 
processes contributing to dissipation in turbulent collisionless plasmas are still an open issue 


and in the recent years it has been argued (see, e.g., Alexandrova et ah, 2013) that below 
the proton scales another turbulent cascade may take place that is described by a different 
power law. It is followed by an exponential law which could be indicative of dissipation. 


In practice, the analysis of in-situ time series provides PSD as a function of the 
frequency in the spacecraft reference frame, P{fsat), that would correspond to Doppler 
shifted wavevector spectra, P{k), under the assumption that the plasma flows over 
the spacecraft much faster than the characteristic time evolution of the nonlinearly 
interacting turbulent spatial structures/eddies (the Taylor hypothesis). In the solar 
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wind the transition between the driving and the inertial range is in general observed at 
freqnencies between 10“^ to 10“^ Hz that wonld correspond to spatial scales related to 
the solar wind correlation/integral length (A) or the typical size of the “energy containing 
eddies” ( Batcheloi] 1970 Matthaeus et ah 1994). The high frequency limit of the inertial 
range in the solar wind and the transition to kinetic regime is marked by a break in the 
spacecraft-frame frequency representation, generally in the vicinity of spatial scales (under 
the Taylor hypothesis) of the order of the proton inertial length or the proton Larmor 
radius (e.g. [Chen et al., 2014). Recently, the variation of this break with the heliocentric 


distance was discussed by [Bruno fc Trenchi] |2014[ High resolution data seem to suggest 
that dissipation may effectively start at higher frequencies, corresponding to the electron 
Larmor radius ( [Alexandrova et ah 2009). The solar wind is a supersonic and super-Alfvenic 
tenuous stream of collisionless plasma emerging from the dynamic solar corona therefore 
discerning “pure” turbulence features from other structures convected from the Sun is still 


an issue (see, for instance, Tu & Marsch, 1995 Bruno et al., 2007 Borovsky, 2008). 

Solar wind observations in the inner heliosphere (between 0.3 and 0.9 AU) suggest 
that the ordering parameter of turbulent properties is the “age” of turbulence, evaluated 
as the time it takes for the solar wind to travel from the Sun to the spacecraft, rather 
than the radial distance. The “aging” of the solar wind turbulence is also characterized by 
a progressive spectral dominance of the “2D” mode of turbulence (characterized mainly 
by perpendicular wavevectors), over the “slab” mode (dominated mainly by parallel 
wavevectors) (Ruiz et al. 2011| . Analyses based on a global mean magnetic held estimation 
concluded that the slow wind exhibits in general features of “2D” turbulence while the 
turbulence in the fast wind is more of the “slab” type ([Dasso et al.j 2005, Weygand et al. 


2011) and the anisotropic state found near the Sun evolves towards a more isotropic state 


at 1 AU. On the other hand, approaches exploring anisotropy through the scale-dependent 
local mean magnetic held (e.g. Horbury et al.| 2008, Podesta, 2009, Wicks et al., 2010 
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Forman et al., 2011) indicate that the high-speed solar wind power spectrum is dominated 


by perpendicular “2D” fluctuations. Simulation results (Chen et al. 2011) show the same 


discrepancy between global and local mean magnetic held approaches and the authors 
conclude that the global mean magnetic held scaling is not able to properly discriminate 


between parallel and perpendicular fluctuations. Smith (2003) shows that at high latitudes 


roughly equal proportions of slab (ID) and 2D coexist in the same plasma element. 

At solar minimum the solar wind is characterized by an increased recurrence of the 
high speed streams (up to 800 km/s and more), lower density and higher temperature, 
whose origin is the meridional extensions of the polar coronal holes. The properties of fast 
and slow wind turbulence were investigated for different phases of the solar cycle, from 


data recorded by Helios (Bavassano & Bruno 1989, 1991 Ruiz et al. 2011), ACE (Borovsky 


2012b), Ulysses (Yordanova et al 2009), Cluster and THEMIS (Weygand et al. 2011), 


Messenger and Wind (Bruno & Trenchi 2014 Bruno et al 2014b). Observations of the solar 


wind by Ulysses at larger radial distances (between 1.5 AU and 5.4 AU), outside the ecliptic 
and close to the solar minimum show that the spectral index of the fast wind inertial range 
turbulence takes values in the range —1.79 < a < —1.55 for the magnetic held components, 
and between —1.52 < a < —1.25 for the total held, \B\. The spectra of the slow wind 
exhibit similar power law behavior but with steeper slopes, —1.95 < a < —1.45 for the 


components of the magnetic held and —1.78 < a < —1.55 for the total held (Yordanova et 


al 2009). It is not clear whether the inertial range spans roughly the same frequency range 


for the fast and slow wind. 

In-situ observations at 1 AU indicate that the median of the magnetic spectral index in 
the inertial range depends on the type of wind: it is shallower for the fast wind (y^ind > 450 


km/s), a = —1.54, compared to the slow wind, a = —1.70, as shown by Borovsky (2012a) 


from ten years of ACE data (1998-2008). The spectral index may take “extreme” values. 
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larger than —1.33 and smaller than —1.95. The steeper spectral slopes are observed at 
1 AU when the solar wind density is larger, the temperatnre is smaller, the speed takes 


smaller valnes and the nnmber of strong directional discontinnities is rednced (Borovsky 


2012a). Close to the high freqnency limit of the inertial range, at proton scales, in the 


vicinity of the flnid/kinetic spectral break, WIND (at 0.99 AU) and Messenger (at 0.42 AU) 
data show that the spectral slope may depend on the power density in the inertial range: 


steeper slopes are observed for larger power in the inertial snbrange (Brnno et al 2014b) 


This freqnency break moves towards smaller freqnencies as the radial distance increases 


(Brnno & Trenchi 2014). 


2. Spectral properties of fast and slow solar wind at 0.72 AU 

We analyze data recorded by Venns Express (VEX) in the solar wind, in the vicinity 
of Venns, at 0.72 AU, between Jannary 2007 and December 2009, dnring the minimnm of 
the solar cycle 24. Since Venns has no intrinsic magnetic field its indnced magnetosphere is 
conhned to shorter distances from the planet, thns VEX spends in the solar wind more than 
20 hours each day of the year. Thus VEX is a unique solar wind monitor that investigates 
the inner heliospheric solar wind on a day to day basis for almost one solar cycle, since 2006. 
The turbulent fluctuations of the solar wind magnetic held considered in this study are 


provided by Venus Express Magnetometer, VEX-MAG, (Zhang et al. 2006)) with a cadence 


of 1 Hz. The data are obtained through an offline calibration procedure by donwsampling 
the 32 Hz resolution data. 

The plasma state (electron and ion spectra and their moments, e.g. density, 
temperature, velocity) is investigated by the Analyzer of Space Plasma and Energetic 


Atoms (ASPERA, Barabash et ah, 2007). ASPERA operates in the solar wind for short 


time intervals of the order of one to two hours, close to the orbit apogee. The ion and 
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electron spectra and their moments are provided with a time resolntion of 196 seconds 
which does not allow a spectral analysis of the solar wind plasma parameters. Nevertheless, 
the estimation of the moments of the ion velocity distribntion fnnction provides the data 
needed to select high and low speed solar wind. Inspired by previous studies, we select 
the slow and fast wind intervals based on a threshold speed value equal to 450 km/s. The 
magnetic held experiment is operating continuously, however we consider time intervals of 
roughly four hours length, close to the VEX apogee and that include the time periods when 
ASPERA is also operating. 
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Fig. 1.— Summary plot of the power spectral density spectrum of Bx, By and B^ from 
VEX magnetic held data recorded between 2007-2009, during the minimum of the 24th solar 
cycle. Yellow and red lines correspond to fast solar wind, gray and blue lines correspond to 
slow solar wind. The black dotted lines mark the frequency interval for which the spectral 
index of each spectrum is computed. An average spectrum is derived as an ensemble average 
of all the spectra obtained for fast wind (the embossed red line) and respectively the slow 
wind (the embossed blue line). The spectral slope of the average spectra is computed for 
the frequency range marked by yellow and magenta respectively. 

Additional constraints on the magnetic held data analysis are imposed by several 






















quality checks. We disregard the time intervals shorter than one hour and that have data 
gaps whose largest size exceeds 30 consecutive points. In general, the total missing data 
points amount to less than 3% of the total number of samples for the selected time intervals 
and they appear to be randomly distributed over time. Thus, the number of data gaps 
in the selected time intervals is small and their length is short. A linear interpolation is 
applied prior to the spectral analysis. From a total of 1094 orbits between January 2007 
and December 2009, only 204 time intervals fulhll the data quality requirements, of which 
48 time intervals correspond to fast solar wind observations {Vwind > 450 km/s). 


The magnetic field components are provided in the Venus Solar Orbital (VSO) 
rectangular frame, with the Ox axis aligned in the Sunward direction and the Oz axis 
perpendicular to the ecliptic plane, in the Northward direction. The PSDs are computed 
with a Welch algorithm (Welch 1967[ ) of averaging periodograms applied on Bx, By, B^, 
and the total held, \B\ for all the selected time intervals. We also compute the trace of the 
spectral matrix of the huctuations. A summary plot of all the power spectra is given in 
Fig.|Tl different colors illustrate different types of wind (fast/slow) and different years. An 
average spectrum is computed as an ensemble average of all the spectra, for the fast and 
the slow solar wind, respectively. 


The spectral power of fast solar wind is systematically larger than for slow wind as 


shown by previous results in other locations of the heliosphere (see the review by Bruno & 


Carbone, 2013 and recent results at 0.38 AU and 1 AU by Bruno et alf 2014b). The power 


spectral densities exhibit a power law regime in the frequency range [5 x 10“^, 10“^] Hz. A 
change in the spectral slope is observed around 0.2-0.3 Hz, close to the local Doppler shifted 
proton gyrofrequency, followed by a frequency range showing a steepening of the spectra 
that could correspond to the lowest frequencies of the non-MHD kinetic turbulent cascade. 
In order to minimize the effect of inherent uneven sampling of the frequency range by the 
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periodogram technique, the spectral index/slope has been computed from PSDs rebinned 
and averaged over equal logarithmic bins of frequency. In addition to the rebinning of 
the frequency range we have also searched for the existence of the power law behavior 
by varying the limits of the htting frequency range. Thus for each of the power spectra 
the spectral index (or slope) is calculated by a linear least-square htting over the interval 
[5 X 10“^, 10“^] Hz. This range would correspond to the high-frequency part of the inertial 
subrange at 0.72 AU and solar minimum. The amount of data intervals is sufficiently large 
to allow a statistical analysis of the distribution of the spectral indices obtained at solar 
minimum, between 2007 - 2009, as illustrated by Fig. 

The a priori separation of fast from slow solar wind allows the independent and 
simultaneous tracking of the evolution of the spectral properties of turbulence for the two 
types of wind. The least-squares linear htting of the distributions of the spectral slopes 
(upper panels of Fig. yields an experimental evidence for the fact that the spectral 
indices do not exhibit any temporal trends over the years while the solar minimum deepens, 
both for slow and fast solar wind. This can be considered an indication that the processes 
that contribute to the power law scaling are mainly local, possibly related to non-linear 
interactions leading to turbulent transfer of the energy between scales. The histograms of 
the spectral indices show that the slopes of the magnetic held power spectra in the inertial 
range ht a Gaussian distribution, as seen in the lower row of plots in Fig. The three 
components of the magnetic held exhibit diherent average spectral slopes for fast and slow 
wind: = —1.57 ± 0.02, = —1.58 ± 0.02, = —1.60 ± 0.02 and respectively 

^slow ^ ^ ^sjou, ^ _^_g4 ^ ^sjow ^ _^_g4 ^ -fast/slow 

the notation for the mean spectral index of the component i and fast or slow type of 
wind. The mean spectral index of the trace of the spectral matrix (not shown) varies from 
_/ast _ _]^_gg 4 _ g Q]^ wind to = —1.65 ± 0.01 in the slow wind. 
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Fig. 2.— The distribution of spectral index for Bx, By and B^ from data recorded between 
2007-2009. Upper raw: evolution in time of the spectral indices computed in the frequency 
interval indicated in Fig. c is the constant (intercept) of the hrst degree polynomial linear 
£t and s is the slope of the fit. Lower raw: the histogram of the spectral indices approaches 
a Gaussian distribution with a mean close to —5/3 for the slow wind and —1.6 for the fast 
wind. 
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Fig. 3.— The spectral index as a function of solar wind speed; we show averages over velocity 
bins of 50 km/s width. Different colors ilustrate different magnetic field components. Data 
is collected in the solar wind by VEX between 2007-2009. The dotted vertical line represents 
the threshold we chose to select fast and slow solar wind. 
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The variation of the spectral index as a fnnction of the solar wind speed is shown in 
Fig. 1^ The average spectral indices of the magnetic held become shallower with increasing 


plasma velocity which is in good agreement with previons work (Chen et al. 2013). The 


continnous trend also seems to indicate that the solar wind velocity may be a controlling 
parameter for the spectral slope. Nevertheless, recent stndies indicate that cross-helicity is 


eqnally important in controlling the spectral behavior (Chen et al. 2013) 


The Venus Express results show similarities with the spectral analysis of WIND 
magnetic held data at 1 AU, during roughly the same time interval, between June 2004 
and April 2009. Data from WIND showed a mean spectral slope between —1.6 (for wind 


speeds larger than 600 km/s) to —1.72 (for speeds smaller than 400 km/s) (Chen et al. 


2013). Another analysis of solar wind turbulence at 1 AU, based on ACE data recorded 


between 1998-2008, suggests a median spectral slope of the magnetic held equal to —1.54 
for speeds larger than 550 km/s and respectively —1.70 for speeds smaller than 450 km/s 


(Borovsky 2012a). Nevertheless, the Venus Express data suggest that in average the slope 


of the slow wind spectra are shallower at 0.72 AU than at 1 AU, while the fast wind shows 
in average comparable slopes. We note the diherences between the mean spectral slopes of 
the three magnetic held components for slow wind conditions and roughly the same mean 
slope for fast wind. On the other hand more spectral power is observed for By and 
components in fast wind. This is in agreement with expectations of hnding more power 
on the perpendicular components than the parallel one as fast wind is more Alfvenic and 
Alfvenic huctuations are mainly perpendicular to the local mean magnetic held and are 
not compressive. At 0.72 AU the Parker spiral (mean held) for fast (700km/s) solar wind 
would be around 24°, not far from the X direction in the VSO reference system. Thus, 


more power should be expected on the Y and Z components of the magnetic held (Klein et 


al. 1993). The anisotropy will be the subject of a future study on the same dataset. 
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3. Discussion and Conclusions 


We have investigated the spectral behavior of the solar wind magnetic held at 0.72 AU, 
for low solar activity, between 2007 and 2009, nsing data provided by Venus Express. The 
Power Spectral Densities of magnetic held components and the trace of the spectral matrix 
indicate that an inertial range of turbulence can be identihed as a power-law behavior, in 
the fast and slow solar wind. More power is contained in the fast wind spectra that also 
exhibit shallower slopes than the slow wind. The mean value of the slope takes values 
around —1.6 ± 0.01 for the fast wind and respectively —1.65 ± 0.01 for the slow wind. Our 
results fully agree with general predictions found in literature about diherent spectral slopes 
of magnetic held huctuations depending on solar wind conditions. 


In particular, (Chen et ah 2013) clearly showed how the spectral index of magnetic 
held and velocity huctuations depend on their Alfvenicity which can be expressed by the 
normalized cross-helicity ac- These authors, analysing 5 years of WIND data, found that 
as ac decreases, magnetic energy starts to dominate on kinetic energy and the magnetic 


held spectrum becomes steeper than velocity spectrum, in line with predictions (Muller 


& Grappin 2004), among others. The tendency for the magnetic spectrum to dominate 


over the kinetic one is a natural outcome for a stationary critically balanced MHD 


turbulence (Goldreich & Sridhar 1995) generated by non-linear interacting Alfven 


waves (Gogoberidze et ah 2012). As a consequence, since we study fast and slow wind. 


which diher in Alfvenicity (Bruno & Carbone 2013), fast wind being more Alfvenic than 
slow wind, we should expect to hnd magnetic held spectra in fast wind less steep than 
those in slow wind. In addition, a comparison with similar data at 1 AU suggests that in 
average the spectra steepen while the slow solar wind is transported between 0.72 and 1 
AU suggesting that non-linear interactions are at work. On the other hand, the fast wind 
exhibits less clear evidence of the radial steepening. 
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The solar wind magnetic field spectral indices between 2007 and 2009 have a normal 
distribution. We do not find a significant temporal trend of the slopes. The average spectral 
slopes of the three magnetic components suggest an anisotropic repartition of power. A 
change in the spectral slope is evidenced in the vicinity of the proton cyclotron radius, 
possibly associated to the transition from the inertial to the kinetic subrange. 


Our results suggest that at 0.72 AU and solar minimum the slow wind exhibits 
in average a spectral behavior closer to the with some differences between 

the mean slopes of the magnetic field components. Thus there are indications that the 
turbulence is anisotropic and models based on the isotropy hypothesis, like Kolmogorov 
and/or Iroshnickov-Kraichnan, are not applicable. Modern theories of anisotropic strong 


MHD turbulence (Goldreich & Sridhar 1995) predict that the perpendicular spectrum of 
turbulence may approach in some conditions the “y-5/3)) gQ^iing. in the absence of resolute 
plasma measurements from VEX we can only suggest that the slow wind spectra are 
perhaps dominated by the perpendicular component. Nevertheless the younger turbulence 
carried by the fast wind is described by shallower spectral slopes that show a tendency 
to approach asymptotically a “y-3/2” power law. This could possibly signify that the 
structure of turbulence exhibits features consistent with models of anisotropic turbulence 


like in Boldyrev ( 2006). In other words, our results may suggest that the slow wind at 0.72 
AU and solar minimum is dominated by filament-like structures at the smallest scales (as 


suggested by Goldreich & Sridhar, 1995), while the fast wind turbulence is dominated by 


the sheet-like structures, at the smallest scales (as suggested by Boldyrev, 2006), possibly 
related to phenomenological turbulent features of the solar wind at the origin, in the corona. 
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